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j^ . Abstract 

Charged current charm and D meson production is studied in detail as a means of de- 
termining the unpolarized and polarized strange sea densities at HERA. All analyses are 
performed in next-to-leading order QCD, including a calculation of the so far unknown 
spin-dependent MS coefficient functions up to 0{as)- It is shown that a decent measure- 
ment is possible in the unpolarized case, provided a sufficient luminosity can be reached in 
the future, while for longitudinally polarized beams it appears to be extremely challenging 
due to limitations imposed on the expected statistical accuracy by the charm detection 
efficiency. 



1 Introduction 

At experimentally relevant Q^ values the different unpolarized ffavor sea quark distri- 
butions q{x,Q'^), q = u, d, and s, are quite distinct, and only for asymptotically large 
values of Q^ they eventually evolve to a common x shape due to the dominance oi g ^ qq 
transitions. These light quark sea distributions are usually treated as massless partons 
in all sets of parton densities, hence requiring some non-perturbative input for their Q^ 
evolutions which has to be determined experimentally. Heavy quarks {rrig ^ ^qcd, i-e., 
q = c and b), however, can be dealt with purely perturbatively, within different methods 
though, which completely determines their x and Q^ dependence, with the heavy quark 
masses rriq being the only physical parameters. 

In recent sets of unpolarized parton distributions it is either assumed that s (= s) has 
the same x shape as u+rf |l|, |^, or s is solely generated by QCD dynamics from a vanishing 
input distribution at some low scale P] , in both cases leading to an overall suppression of 
the X integrated second moment J dx xs in the light sea of about 50% at Q^ — 5 — 10 GeV^ 
[4-7], presumably due to the larger mass of strange quarks. The entire x dependence of the 
flavor decomposed unpolarized light sea is, however, still rather uncertain. While some 
information about u and d is now available from various source^, and all data indicate 
that d is greater than u, s can be inferred only from CCFR data on deep inelastic neutrino 
charm production [^ ^ for the time being, with a NuTeV update to be expected in the 
near future p. An alternative extraction along ^F2^{x,Q'^) — SF^* {x,Q'^) ^ xs{x,Q'^) 
(or equivalently from corrections to the F^ / F!^^ ~ 5/18 rule) combining CCFR and 
NMC [|lOl data cannot be reliably performed because it suffers from the fact that s emerges 
only as a small residual of two large numbers {xs <^ F^ , FI^^)^ which furthermore appear 
to be incompatible at low x [|ll|. At present, results for 1^2^'^ — SF^* seem to be in conflict 



with the CCFR charm production data [11, 12], and if this tendency persists with future 



NuTeV data, it requires further clarification [^ 



The leading order (LO) contribution to charged current (CC) charm production in 



^Recent compilations can be found, for instance, in Refs. 



deep inelastic scattering (DIS) is given by the 0(0°) parton model process 

ly+s'^c , (1) 

depicted in Fig. 1(a), where s' denotes the Cabibbo-Kobayashi-Maskawa (CKM) 'rotated' 
combination 

s' = \Vcs\^ s + \Vcdf d (2) 

with \Vcs\ = 0.9745 and \Vcd\ = 0.2205 [0. Due to the smallness of \Vcd\^ in (D the 
process (^ is expected to be essentially sensitive to the strange sea content. Only at large 
X, where quark sea contributions are less relevant, the \Vcd\ suppression is balanced by 
the valence enhancement of the well-known d^ix). 

In next-to-leading order (NLO) QCD, however, this simple picture is spoiled, and the 
complete set of Feynman diagrams in Fig. 1 has to be considered. Apart from the virtual 
and real 0{as) corrections to ([I|), the genuine NLO gluon induced subprocess 

W+g -^ cs' (3) 



has to be taken into account as well, which may yield a significant contribution [|14 
to the charm production cross section, hence representing an important 'background' 
for any extraction of the strange sea. In case of inclusive charm production these NLO 
corrections have been calculated for unpolarized target nucleons in different regularization 



prescriptions such as the conventional MS scheme [|^, [T2[, which we henceforth adopt, 
or the ACOT scheme [[161 , P^ , which also takes into account possible effects of a finite 
strange quark mass. Recently these fully inclusive calculations were extended to the 
experimentally relevant case of momentum (z) distributions of the produced D mesons 
[18-20]. Such detailed production cross section considerations seem to prefer a softer, 
'heavy quark-like', strange sea [§[ over s oc m + Jinputs P, |l], H, but further experimental 
clarification is certainly highly desirable. 

Besides new fixed target neutrino data from NuTeV [Q, one interesting possibility to 
shed more light on the strange density and its Q^ evolution would be, of course, to study 
CC charm, i.e., dominantly D meson, production in e^p collisions at HERA, provided 



a sufficient luminosity can be reached. We shall perform a closer analysis of this option 
and the impact of the gluonic 'background' (0) in this case below. It should be noted 
in passing that a high precision measurement of CC charm production in e~p and e~^p 
collisions could possibly reveal also the relevance of recent claims |2l| that strange and 



anti-strange densities may differ, i.e., s y^ s, contrary to what is assumed in all analyses 
of parton densities so far. 

Turning to longitudinally polarized parton densities A/, defined by 

Af{x,Q')^U{x,Q')-f4x,Q') , (4) 

where /+ (/_) denotes the distribution of a parton / with its spin (anti-) aligned to the 
parent nucleon's spinQ, much less is known experimentally about the flavor decomposition 
of the light sea or even the gluon density Ag. Information on the A/ is so far almost 
exclusively available from fully inclusive polarized DIS, i.e., structure function measure- 
ments [^ , which are only sensitive to specific non-singlet and singlet combinations of the 
spin-dependent quark densities and not to a full flavor separation or Ag. Thus all current 



sets of polarized parton densities, such as, e.g., the GRSV |Q or GS ||2J] distributions, 
have to fully rely on certain assumptions when providing flavor decomposed quark den- 
sities, which are often biased by unpolarized measurements and, of course, remain to be 
checked. 

The knowledge of the flavor decomposed polarized densities, more specifically of their 
ffist moments Af{Q^) (as obtained by taking the x integral in (^)), is moreover required 
to understand how the nucleon's spin S^ = 1/2 is shared by its 'constituents' as a function 
of the momentum transfer Q^, 

S^ = l = l E [^^(^') + ^^"(^')] + ^9{Q') + LUQ') + mQ') , (5) 

q=u,d,s 

where L| (Lf ) denotes the orbital angular momentum contribution of the quarks (gluons)^. 

^By taking the sum instead of the difference in (||) one recovers the unpolarized (hehcity-averaged) 
parton densities /. 

•^Of course, in NLO the decomposition on the right-hand side (r.h.s.) of ^ becomes factorization 
scheme dependent, and one always has to specify the scheme one is referring to when quoting values for 
the first moments Af{Q^) or Lf^. 



The presently available semi-inclusive spin-dependent DIS measurements [^ are still 
not conclusive enough to disentangle different flavors reliably, but some progress has to 
be expected soon in particular from the HERMES experiment. Together with upcoming 
measurements of W boson production at the polarized BNL-RHIC pp collider this may 
yield some information about An and Ad. However, a direct measurement of As like in the 
CCFR neutrino DIS experiment turns out to be extremely remote despite that neutrinos 
have deflnite helicity, since tons of target material would have to be polarized. Thus other 
possibilities have to be examined here. Since it is no longer inconceivable that HERA can 
be operated at some stage in the future in a longitudinally polarized coUider mode PB| , 
it was suggested to study CC charm production to decipher As |^ along similar lines as 
discussed above for unpolarized e^p collisions. However, these studies neither have been 
performed in NLO, nor do they include a realistic estimate of the expected statistical or 
theoretical uncertainties for such a measurement at a polarized HERA. It is the main 
purpose of this paper to provide the complete NLO framework for CC inclusive charm 
and momentum {z) differential D meson production in the MS schemeQ with polarized 
beams and to study the prospects of a measurement of As at a polarized HERA. 

The remainder of the paper is organized as follows: in Section 2 we outline all relevant 
technical details of the calculation of unpolarized and polarized CC charm production 
in NLO, mainly focusing on additional complications which arise in the momentum (z) 
differential and in the spin-dependent case. Section 3 is devoted to a detailed numerical 
analysis. First we discuss the prospects of a determination of the unpolarized strange 
density at HERA, then we turn to the polarized HERA option and its potential of learning 
something about As. The Appendix contains the polarized coefficient functions, which 
are our main analytical results and are too long to be presented in the text. 



"'We compare our calculation to existing inclusive MS |Q as well as fc™'"- [g^ and mass-regulated 
1 results at the end of Section 2. 



2 Technical Framework 

To derive the cross sections for inclusive charm and momentum z differential D meson 
production in longitudinally polarized CC DIS we follow closely the corresponding unpo- 
larized NLO calculations in [O, O, Il8[. Since only very few technical details have been 



presented in [^ in the experimentally more relevant case of heavy meson production, we 
shall give a brief overview of the most important, non-trivial calculational steps as well. 
We mainly focus, however, on the complications arising due to the appearance of 75 and 
the Levi-Civita tensor e^j^^p^, in course of the calculations, which requires special attention 
in the polarized case as we will discuss below. 

The NLO corrections to the LO parton model CC production mechanism^ ([l|) stem 
from the boson gluon fusion (BGF) and real gluon emission subprocesses (^ and W^ s -^ 
eg, respectively. In the latter case also virtual corrections to (|l]) have to be included. 
All contributions are represented by their Feynman diagrams in Fig. 1. As usual, the 
NLO {0{as)) diagrams comprise soft and coUinear divergences, and one has to choose a 
consistent method of regularizing these singularities so that they become manifest. For 
this purpose we work in the well established framework of dimensional regularization in 
n = 4 + 2e spacetime dimensions. The divergences then occur as poles ~ 1/e and ~ 1/e^ 
in the physical limit n — i> 4. The latter double pole terms only arise in the quark initiated 
subprocess when soft and mass/collinear singularities coincide. While these double and 
the single poles from soft gluons in virtual loops and real soft gluon emission have to 
cancel by the KLN theorem, there remain mass/collinear poles ~ 1/e in the gluon and 
quark initiated NLO corrections stemming from the region in phase space where a strange 
quark propagator goes on-shell. This can happen either when the initial state gluon splits 
into a collinear ss pair or when the initial state strange quark radiates off a collinear 
real gluon. These poles have to be removed from the production dynamics by factorizing 
them off into the renormalized (scale-dependent) strange sea density, where we adopt 
the commonly used MS factorization prescription. Charm quark propagators in collinear 



^ For simplicity we ignore the CKM s <-> d mixing (0) in this technical section. It is, however, properly 
taken into account in all phenomenological applications in the next section. 



g -^ cc, c —>■ eg subdiagrains are protected from going on-shell by the (heavy) charm 
quark mass, which thereby acts as an effective cut-off of non-perturbative long distance 
strong interactions. 

In the following we give a unified description of the calculation of the relevant produc- 
tion processes for unpolarized and longitudinally polarized initial states. The possibility 
to obtain the unpolarized results 'simultaneously' provides a useful check of the correct- 
ness of our results, and we fully agree with the unpolarized results presented in [|^, [1^, |1^ . 



To be specific, we calculate the contributions of incoming quarks and gluons to the un- 
polarized and polarized structure functions H!^'^ and AHf'^, respectively, as depicted in 
Fig. 1 by properly projecting the hehcity dependent squared matrix elements |A^^'^(±)p 
onto the structure function i: 

( TT<i,9 ^ ( p'^f^ 

I i 



\AHr 




^^•^(+)r±i-M^i(-)r dPS2 . (6) 



The indices a and /? in (||) indicate the polarization indices of the W^ boson, and dPS2 
is the two body phase-space as defined in (0) below. The operators P" are given in 
Eqs. (B9)-(B11) of [|T5| and project for i = 1, 2, and 3 onto the relevant unpolarized 
structure functions Fi, F2, and F3, respectively. The structure functions F4 and F5 in 



15| do not contribute to the lepton-nucleon CC cross section if one assumes a vanishing 
mass for the lepton as we will do in the following. Since the polarized structure functions 
appear in a similar way in the hadronic tensor as the unpolarized ones (cf. also Eqs. ( ]T9| ) 
and pop below), the same projection operators apply in this case if one identifies 

APf = Pf , APf = Pf , APf = Pf . (7) 

The operators AP" in ([^) then project onto the relevant polarized structure functions 
gi, g3, and 5^4, respectively, and other structure functions like g^ and gj again do not 
contribute for a vanishing lepton mass [^. The projection onto the helicity states h = ± 
of the incoming strange quark or gluon in the matrix elements A^ in (|^) is achieved by 
using the standard relations (see, e.g.. 



u{ps, h)u{ps, h) = - ^sil - h'j5) 



for incoming strange quarks with momentum ps (analogously for anti-strange quarks) and 



^f^iP9,h)el{pg,h) = - 



Pn^ 






(9) 



Pa-<1. 
for incoming gluons with momentum Pg [q denotes the four-momentum of the W^ boson). 

The presence of 75 and the totally anti-symmetric tensor e^^p^ in the V-A vertices 
and P3 = AP" , respectively, and - in the polarized calculation - also due to (H) and 
(y), introduces some extra complications, because their purely four- dimensional origin 
allows for no straightforward continuation to n 7^ 4 dimensions. We choose to handle 
these quantities in the HVBM scheme [^, which was shown to provide an internally 
consistent continuation of 75 and e^^ypo- to arbitrary dimensions. This prescription is also 
implemented in the package TRACER |33], which we use for all trace calculations in n 



dimensions. In the HVBM scheme ||32[ the e-tensor continues to be a genuinely four- 
dimensional object, and 75 is defined as in four dimensions, implying {7^,75} = for 
/i = 0, 1, 2, 3 and [7'^, 75] = otherwise. This effectively splits the n dimensional space 
into two subspaces: one containing the four space-time dimensions and one containing 
the remaining {n — 4) dimensions, denoted as 'hat-space' henceforth. The price to pay 
is that in the matrix elements squared in (^ we then encounter not only conventional 
n dimensional scalar products of two momenta, which can be expressed in terms of the 
usual partonic Mandelstam variables, 

s = {Ps,g + qf , t = {q- pcf , u = {ps,g - Pc)^ (10) 

in our case, but also similar scalar products in the hat-space. However, in the parton-boson 
cm. system with the incoming momenta chosen to lie in the ±z direction, all possible 
(n — 4) dimensional scalar products of the two final state momenta can be expressed in 
terms of a single hat momenta combination k^ = —k ■ k due to momentum conservation. 
The k"^ terms do not contribute to the unpolarized calculations^ |T5|, |^, |18|, but are 



important for polarized DIS due to the additional appearance of 75 and e^i^p^ in Eqs. 
and (|9|) as we will discuss in more detail next. 



^In principle these calculations can be performed also by using the theoretically inconsistent anti- 
commuting 75 prescription of p4[ in n dimension as was done in Pq|. 



The partonic two particle phase space dPS2 for one massive and one massless parton 
in (0) is given by 



(3-m2)2 



^PS2 = i- (47r)-^^^ Vn dy dp (e) (11) 

Stt s T{e) Jo Jo V / 

1 s-m? 1 r/„ ^2N2n^ .i 



i7T s Tl+e) 



m2)2 



, ,.v- ... -. - vl2) 

In (|T2|) the integration over d/c^ has been carried out, and the standard n dimensional 
phase space [|1^] is recovered. This can only be done if either the matrix element squared 



\A4\ in (^ trivially does not depend on k"^ or if the k'^ dependent terms in |A1| are not 
multiplied by a sufficiently singular term ~ (1 — y)~^, because the extra subvolume of the 
phase space, which is available for k"^, is intrinsically of order 0{e) due to the l/r(£:) in 
(p!T|). These conditions are always met in the unpolarized case. 

Eq. ([TI|) leaves the two options to either fully integrate the {A)H^'^ in (|^) over the 
entire phase space or to consider more differential observables which can be obtained 
from some Jacobian according to d/dS = d/dy dy/dE, where S stands for any kinematical 
observable that can be expressed hjy= (l+cos6'*)/2, with 9* being the VT^-parton c.m.s. 
scattering angle. A more detailed discussion including subtleties arising from endpoint 
(soft) singularities will be given below when we consider the fragmentation spectrum of 
charm quarks in CC DIS. 

I I 2 

The helicity-dependent matrix elements A1^'£(±) in (|^) can be easily derived from 
standard Feynman rules and will not be given hereQ. They can be conveniently expressed 
in terms of the partonic Mandelstam variables (0), which in turn can be written as 

t = zl^s + Q')is-ml)il-y) (13) 



s 

2 

2 
c ' 



u = [ml + y{s — ml)~\ + m 



where Q^ = —q^ and 

2ps,g • q V Q J s + Q^ 

^The non-trivial virtual corrections are explicitly calculated in Ref. |l3] which we confirm after elimi- 
nating a misprint in the coeflicient A2 in Tab. 1 of |l5|, which should read p^ Ka instead of Ka/2. 
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is the partonic analogue of the slow rescaling parameter ^ = x{l + ml/Q'^) |^. Its 
introduction is required in NLO for a consistent factorization of initial state collinear 
singularities [l^ . 



Within dimensional regularization the soft and collinear singularities in {A)Hf'^ can 
be isolated using standard distribution- valued expansions |^ of the type 



r(i-y) 



-l+e 



1 



S{l-y) 



1 



(1 - y)- 



+ e 



In(l-y) 

i-y 



+ 



i-y 



(15) 



where the '+ '-distributions are defined in (|A4|) in the Appendix. The relevant expansions 



for the initial state variable ^' can be found in Eqs. (30)-(32) of |T^. After isolating 
soft and collinear poles and canceling the soft poles we can - according to our discussion 
below Eq. ([TID - transform the results to the final state charm quark momentum scaling 
variable ( =Pc- Ps,g/q ■ Ps,g using 

dy s 



(16) 



s — m^ 

Care has to be taken when applying the Jacobian in Eq. (|T^) to distribution- valued expres- 
sions in the variable y as in Eq. (^5|). There the transformation does not merely amount 
to a multiplication with dy/dC, but rather has to be found by changing the integration 
variable when folding the distribution with some test function such that: 

5{l~y)^5{l-0 ; TT^^TT^ ; "T r - ^r^ ,(17) 



(l-y)+ (l-C)e 

where the '©'-distributions are defined again in ([A4|) , A 
(1 — A)^'/(l — A^'). Furthermore, because Cmin — * 1 as ^' - 



rin(l-y)l 

[ 1-^ J 


+ 


rln{l-C)l 


-'- Simin 


_i-C_ 



g2/(g2 ^ ^2)^ ^^^ ^^ 



5(1 -0/(e', = 5(1 -0^(1-0 



daf{i\a) 



S,min 



«'=i 



where / may be either a function or distribution. We note here that by an analogous 
Jacobian transformation as in Eqs. ([T5|)-(p!8D we could in principle obtain differential 
distributions in the transverse momentum p^ (= a/s/4 — ml sin 6*) of the produced charm 
quark as well. 

Using the standard tensor decomposition of the hadronic tensor, the structure func- 
tions F^^ 2 3 (unpolarized) and g^^ 4 ^ (polarized) refer to the following double [triple] 



differential CC e^p cross sections 



dx dy [dz] 
dx dy [dz] 



G2 a 



27r(l + Q'^/M^ ^2 



'w) 



y)Fr + y'xFr ± y{l - ^)xFr 



y. 

2' 



r^2 c 



27t(1 



QVM^ ^2 



'w) 



y)g. 



w^ 



2/x w^) 



4 +y xg^ ± y{l - -)xg^ 



(19) 



(20) 



for a fully polarized lepton beam (left-handed e~ or right-handed e"*") scattering off an 
unpolarized (|19|) or a polarized ( pO]) target and where Gp, Sep and Mw denote the Fermi 
coupling, available c.m.s. energy squared, and VT-boson mass, respectively. The longitu- 
dinally polarized cross section (iAa in (^) is defined as the difference [da^ — da^), where 
=^ (^) denotes the direction of the proton (lepton) spin. In case the incident lepton is 
not fully polarized, i.e., Pg 7^ 1, the r.h.s. of (p!9| ) and ( pOD have to be multiplied by an 
overall factor (1 + Pe)/2, which amounts to a factor 1/2 for an unpolarized lepton beam 
(Pe = 0). One should note the change of sign of the P3 and gi terms in (p!9|) and (pO]) , 
respectively, when the electron e~ is replaced by a positron e"*". In Eqs. (|TP|) and ( PDj ) x 
and y are the standard kinematical deep inelastic variables (Bjorken scaling variable and 
inelasticity, respectively) and, moreover, z = po ■ P/q ■ P {pd being the charmed hadron's 
momentum) is a final state scaling variable. 



The CC structure functions entering d{A)a^ ^ in Eqs. ( pJl) and ( poD are obtained by 
the following convolutions^: 






Tf 



x,Q^ 



^x, z, t^ 



where 






As ) Zn 



+ 
+ 



+ 



/ 



m 



9 r^ .2 



^^^.(e',.l,A)||;^(|,,^) 



^-(e,/i|)| De{z) 

271" Jg i' Jmax(z,Cmin) C 

Hf 



W 



{^'^,(e',c,/i|,A) 



^^, (r, C, /iF, A) 1^ ^ ^J^,f^l) 



A^^e 



/^.(^) . 



-9i 






1 r-P3^ 

9\ 






aVI'^^) 



(21) 



aVI'^^) 



(22) 



1 
2^ 1-^2 



(23) 



^With obvious adjustments for da'^ ^ in Eqs. ( |2l| ) and ( p^ ) below: (A)s ^ (A)s and {.^3,^/3 4} 
-{^3,^3,4}- 
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The superscript 'c' in (|21|)-(p3|) indicates that we restrict ourselves only to the charm 
production contribution to deep inelastic CC structure functions and the cross sections in 
(|19|) and (|20|) . For simplicity we set the factorization scale ^p equal to the renormalization 
scale /iij in (|2T|) and (p2D and fix both at fij^ = fij^ = Q'^ + ml. The coefficients Hf'^ can 
be found in P^ , [T^ while the polarized AHf^ are new and listed in the Appendix. Please 
note that the Aiff (^, (^,/i|,, A) in the Appendix (exactly as the iff (^, C,/i|', A) in p!8| ) 



comprise terms of the type [/(0]+fl'(0 where g is a. singular function at ^ = 1. These 
terms seem to be ill-defined at first sight. They are, however, completely well-behaved 
on the phase space of the double convolutions in Eq. (|22|) because the integration volume 
A^ = 1 — ^min vanishes at ^ = 1 such that [g{^){l — Cmin)]€=i is always finite. 

In Eq. (^) the (factorization scale independent) charm fragmentation function is taken 
as [|3 

D,{z) = n[z[1- Z-' - ej{l - z)fy' (24) 

with the normalization constant A^ being related to e^ via J„ dzDc{z) = 1. The 'hard- 
ness' parameter will be fixed to be Ec = 0.06 for our numerical applications in the next 
section in agreement [|T^, ^ with fixed target neutrino data [4-6]. For our phenomeno- 



logical considerations in the following section the precise value of e^ is, however, of minor 
importance. 

An important comment about our final expressions for the Aiff given in (|Al|)-( |A3|) 
should be made. A naive calculation, as outlined above, gives: 

A^t3,4,i(^' C, ^l, A) = ^ti,2,3(e, C, /il, A) - 4C^(1 - 05(1 - C) . (25) 

The difference 4Cf{1 — C,)S{1 ~ () stems, however, from a too naive factorization of the 
initial state coUinear s -^ sg singularity in the polarized case in ra 7^ 4 dimensions. In 
order to restore helicity conservation at the strange quark-gluon vertex [Q the finite 
renormalization 

A^t3,4,i(^' C, /il, A) -^ A^t3,4,i(^' C, /^l, A) + 4C^(1 - 0S{1 - C) = i^ti,2,3(e, C, f^l, A) 

(26) 
has to be considered, leading to our final results in the Appendix. These NLO (MS) 

11 



coefficient functions AiJ|4 ^(^, (^,yu|., A) hence coincide with the corresponding unpolar- 
ized expressions if^2 3('C) C)/^F) '^) HH) as must be due to the same tensorial structure at 
the partonic level, the nature of the V-A interactions, and helicity conservation at the 
(massless) strange quark-gluon vertex. It should be remarked that the fully inclusive 
quark coefficients AH^ {C^ , fi^ , X) have been already obtained in [^, and we fully agree 
with these results. 

For the H? the relation between the small strange quark mass limit of the fully massive 
{tjIs^c 7^ 0) BGF process and the MS results of |l^ has been established in Eq. (5) of 
Ref. ||19[. We note here that analogously the AH^{^, (, /i|,. A) in ( |A^ ) can also be obtained 
from the m^ — > limit of the general, fully massive spin-dependent BGF expressions 
in Eq. (10) of |^. This is a non-trivial cross check of our MS results. Similarly our 



expressions for AHf{^, fi'jp, X) in (|A1 1|) agree with the fc™™ (transverse momentum cut- 



off) regulated results in Eq. (13) of [^ for A;™™ = in the limit m^ — > 0. Finally, taking 
the limit rric — *> of our AH^'^{$^, ,/x|,. A) we recover - apart from obvious coUinear logs 
- the massless results in ||3^ after transforming them to the conventional MS scheme 
by setting the Af^'^ = in |39|. The results in |3^ originally refer to a scheme with a 



non- vanishing gluonic contribution to the first moment of gi, g^ and g^ as do the A;™' 
regulated results in Eq. (15) in |^9|| when taking the massless limit. 



3 Determining s and As at HERA 

Equipped with the required technical framework, we now turn to a detailed discussion 
and numerical analysis of the prospects of determining s and As at HERA in the future. 

To begin with. Fig. 2 shows the unpolarized z differential D meson production cross 
section e~p —>■ DX according to Eq. (|19|) for a fixed value of x = 0.05 and where we 
have integrated over Q'^ = SepXy > 500 GeV^ employing the cut 0.01 < y < 0.9, with 
,^^/S^ = 300 GeV for HERA. Unless otherwise stated we use the GRV-94 0] distribu- 



tions in all our unpolarized calculation^^. The different contributions to the NLO cross 



''This is mainly to avoid any inconsistencies due to different Aqcd values, spurious violations of the 

12 



section da /dxdz are shown separately in Fig. 2 to demonstrate the impact of the gen- 
uine NLO gluon induced subprocess. Needless to say, only the total NLO cross section 
is a physically meaningful observable. As can be seen, the gluonic contribution becomes 
increasingly important and eventually dominates at small values of z, where the Born 
approximation becomes completely meaningless. The sharp rise for 2; ^ can be traced 
back to the 1/C behaviour in the g ^ cc splitting in the w-channel subprocess depicted 
on the r.h.s. of Fig. 1(d). This behaviour should not be considered as a destabilization 
of the perturbative expansion by NLO corrections because it is entirely due to the first 
contribution from charm quarks produced by strong interaction dynamics which one even 
should expect to be important at high Q^. Note that in the region of BGF dominance 
around z < 0.1 the production dynamics and the resulting steep shape is very similar to 
the neutral current case [^ , where the gluon fusion production channel is - within fixed 
order perturbation theory - classified as leading order. To avoid any confusion we will, 
nevertheless, in the following continue to count all 0{as) contributions as 'NLO' in the 
CC case under consideration here. Let us note that, on top of an s{x) measurement, the 
mere observation of a rising cross section towards small z at HERA would be an inter- 
esting experimental confirmation of the underlying QCD dynamics and evolution effects 
because such a behaviour is completely absent at fixed target scales [18,4-7], where the 
contribution from the charm quark in the u-channel of the BGF is small. Integrated over 
z, the steep small z rise results in a quasi-coUinear logarithm ~ InQ'^/ml, which may 
lead to substantial gluonic corrections in the inclusive cross section da/dx unless they are 
removed by acceptance losses or suitable cuts, as will be discussed in more detail below. 

Both (gluon- and quark-initiated) NLO corrections are sizeable in the entire z range 
as can be inferred from comparing the Born result with the individual NLO contributions 
in Fig. 2. The NLO corrections tend to soften the D meson z spectrum and shift the 
peak due to the Peterson fragmentation spectrum (0), located at around z ~ 0.75 — 0.8 
in the Born term, towards smaller values of z. However, we have used NLO distributions 



positivity requirement |A/| < /, etc., in the calculation of the polarization asymmetries (^7|) below, 
since the spin-dependent GRSV parton distributions |g^, which we adopt, were being set-up with the 
unpolarized GRV-94 distributions |40| as reference. Anyway, the results obtained using the recent GRV 
update P] are almost indistinguishable. 
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and a NLO value for Sc in (|2^ ) for the calculation of the Born term, and hence a LO and 
NLO comparison does not reflect a real K factor, i.e., the ratio daNLo/d'^LO^ which has 
to be obtained in a consistent LO and NLO calculation using LO and NLO densities and 
parameters, respectively. We refrained from doing so here, since on the one hand possible 
differences between the LO and NLO strange densities are hardly known yet, and on the 
other hand our main purpose is to illustrate the impact of the NLO corrections undiluted 
by different choices of the parton densities and other parameters. Furthermore, K factor 
considerations are somewhat misleading for the differential observable considered here 
because the 'LO' Born term is at the partonic level by construction^ sharply peaked in 
the forward (proton) direction ~ 5(1 — C) oc 6{1 + cos9*), and a continuous spectrum 



in z is only achieved by smearing the delta peak with the Peterson function in (p4D, i.e., 
purely by non-perturbative hadronization effects. Only in NLO a differential distribution 
which covers the full phase space is obtained already at the (perturbative) partonic level 
and then translates into a realistic hadronic D meson momentum spectrum. It should be 
mentioned as well that in the very large z region, z > 0.85, perturbatively large ln(l — z) 
terms in the NLO subprocesses have to be resummed |^, and non-perturbative higher 



twist contributions will perhaps become relevant |^ such that our results in this region 



cannot be trusted (as is obvious, of course, from the negative cross section obtained here). 

In Fig. 3 we compare the results for the z integrated D meson production cross section 
da/dx for two cases: either integrated over the entire z range (inclusive cross section) or 
restricted to the region 0.2 < 2; < 1, which turns out to be a suitable theoretical cut to 
strongly reduce the gluonic 'background' and perhaps also mimics a poorer experimental 



resolution a low z [Q. As expected from our results in Fig. 2 we obtain again large NLO 
corrections in the first case, mainly due to the sharp rise for z ^ in NLO, whereas in 
the second case the corrections cancel to a large extent. This cancellation can be readily 
understood from the results shown in Fig. 2: both NLO contributions change sign at 
around 2; ^ 0.5 — 0.7 and roughly integrate to zero in the chosen z range|^ 0.2 < z < 1. 

^°It sums up the leading collinear logs. 

-'^-'^One can restrict the integrations also to the range 0.2 < z < 0.85 to avoid the region in the cross 
section where ln(l — z) logarithms become dominant. This would lead only to a slightly larger NLO 
result. 
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The small-z region where the gluonic correction becomes dominant is simply left out, 
thereby minimizing the residual NLO effects. It should be mentioned here that in the 
charge conjugated case e~^p -^ DX the d^ ^ c valence enhancement at large x, which we 
have briefly discussed in the Introduction, becomes significant at the ~ 10% level around 
X ~ 0.1, but dominant only for very large x > 0.5. 

A comment should be made, however, about the NLO corrections which seem to be 
anomalously large at small values of x in the fully inclusive case (0 < 2; < 1), where 
they are dominated by a large quasi-collinear logarithm hiQ"^ /m?^. One may attribute, for 
the moment, this logarithm to a corresponding W~c -^ s contribution^ by introducing 
a resummed, massless charm density c{x). The large NLO corrections seem then in 
turn to imply a larger contribution from charm than strange quarks, i.e., c{x) > s{x) in 



such a language, see also Fig. 2 (r.h.s.) in [^, which would be theoretically not very 



appealing and would therefore seriously question the perturbative reliability of the fixed 
order logarithm InQ^/m^. One has to be careful with such interpretations though, since 
for small values of x and hence large values of y of about 0.7 at HERA, the structure 
function F3 becomes important in (|T9|) . In F^^ s and c enter with negative and positive 
signs, respectively. Taking into account the x and y dependent weights in front of Fi, F2, 
and F3 in (|1^,) it turns out that the negative s contribution in F3 cancels to a considerable 
amount the positive contributions in Fi and F2, whereas all charm contributions add up. 
Therefore the surprising result that charm quarks contribute a larger portion to da than 
strange quarks (mimicing 'c(a;) > s(x)') is simply an effect of the electroweak couplings 
in that particular kinematic region. 

Another indication that the large gluonic correction ~ InQ'^/ml in Fig. 3 does not 
necessarily call for coUinear resummations comes from an 0{al) calculation of CC heavy 
quark production in the asymptotic limit Q^ ^ ml |^5[ . Here it was found that although 
the 0{as) gluonic contribution is large, the O^a^) terms hint at a completely stable 
and well-behaved fixed order perturbation series. As an indication of the perturbative 
reliability of our results it should be furthermore noted that the dependence of our results 

^^ Obviously such a process would not be considered as CC charm production since there is no charm 
activity in the final state (see also the discussion below). 
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in Figs. 2 and 3 on the precise value of the factorization scale fip is rather weak. Variations 
of fj,p in the range 0.1{Q^ + ml) < fip < 10{Q^ + m^), i.e., by two orders of magnitude 
around our default value fi'jp = Q"^ + m1^ change the results for da/dxdz and da /dx shown 
in Figs. 2 and 3 by at most ~ ±10%. On top of this comes a more practical reason for our 
preference for fixed order perturbation theory in the CC case. As mentioned above the 
resummation of quasi-collinear logs ~ InQ^/m^ requires simultaneously the introduction 
of a masslessly evolved charm density c{x) entering a W~c -^ s production channel, where 
the corresponding c from the gluon splitting must be thought of as hiding in the hadronic 
debris. This means that we loose any information whether the event will be tagged as 
a charm event experimentally. Even if we assume that the c escapes always undetected, 
we have no possibility for a gauge and renormalization group invariant separation of this 
production channel from the (tagged) rest of the cross section. On the other hand, in a 
differential fixed order calculation we can, as demonstrated above, apply suitable z —^ Q 
(or p^ -^ 0) cuts to exclude the quasi-collinear region. 

It should be mentioned that the x value in Fig. 2 has been chosen to guarantee a 
sizeable cross section, i.e., a realistic chance to actually measure s. For larger values 
of x, da/dxdz shows qualitatively similar features as the ones illustrated in Fig. 2 but 
at a much reduced cross section as can be inferred from Fig. 3. Of course, the gluon 
contribution becomes less important with increasing x, and the rise for 2; — * becomes 
much less pronounced. With the upcoming luminosity upgrade at HERA it should be 
feasible to study these CC reactions with sufficient statistics. The obtained cross sections 
in Figs. 2 and 3 are all in the ball-park of about 30 — 50 pb for x ~ 0.05 and thus seem to 
be measurable, assuming an integrated luminosity of £ = 200 — 500 pb~ in the future, 
even if one takes a rather small charm detection efficiency of about 1 — 2% into account. 
The charm detection is actually the limiting factor for such measurements, in particular in 
the polarized case as will be shown below. The z-integrated cross section including a cut 
z > 0.2 in Fig. 3 is particularly suited for a measurement of the strange density at HERA 
since the 'background' from BGF drops out almost completely, and the full NLO cross 
section can be approximated by its Born term, which is directly sensitive to s. However, 
within the achievable accuracy it seems to be virtually impossible to distinguish between 



16 



different currently available strange densities as is illustrated in Fig. 4. The differences 
between the results obtained using the GRV-94 densities [^] (or the recent update 0) 
or the MRST distributions 0] are not very pronounced at Q^ > 500 GeV^ since they are 
washed out by the Q^ evolution. A possible improvement would be to add the results 
for D and D meson production, i.e., the results obtained for e^p and e^p CC reactions, 
which would obviously double the statistics. The price to pay is of course that one would 



loose any sensitivity to possible differences between s and s ||21[| . However, it seems to be 
difficult anyway to test the latter at HERA unless s and s would drastically differ which 
seems to be not very realistic. To finish this discussion it should be stressed that the 
differences between the GRV and MRST results in Fig. 4 are indeed mainly due to the 
different assumptions about the strange density. Any differences in the gluon distribution 
are not important in that particular kinematic region (in Fig. 4(b) the gluon contribution 
almost cancels anyway), for instance, using the 'large' gluon or 'small' gluon set of MRST, 
MRST(gt) and MRST(gJ,) [|l|, respectively, instead, hardly leads to any changes. 

Let us now turn to the polarized case. As already mentioned in the Introduction, 
such measurements could be performed at HERA as well provided that the option to 
longitudinally polarize both beams |^ will be realized in the future. Since large polarized 
targets as required for neutrino DIS are not likely to be ever build, it should be stressed 
that HERA would be a unique place to study CC DIS (among other reactions which can 



be only analyzed at a polarized ep collider, see ||26||). 



Figs. 5 and 6 show the polarized z differential and integrated D meson production 
cross sections according to Eq. ( ppj ) in a similar way as above in the unpolarized case in 
Figs. 2 and 3, respectively. Unless otherwise stated we use the GRSV 'standard' set of 
spin-dependent parton densities [^ in our calculations, which provides a rather large. 



negatively polarized strange density, i.e., As{x, Q"^) = As{x, Q^) < for all values of a; and 
Q^. For the differential cross section dAa/dxdz in Fig. 5 one observes the same qualitative 
feature as in the unpolarized case: the NLO corrections become increasingly important 
for small values of z due to sharp rise of the gluonic contribution. However, the polarized 
gluon density Ag is only very weakly constrained by presently available data, and hence 
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the actual size of the gluonic correction is extremely uncertain and strongly dependent 
on the chosen set of parton densities. A meaningful measurement of As is therefore only 
possible if /S.g becomes better constrained in the future, which is not unlikely in view of 
the forthcoming experiments at BNL-RHIC and CERN (COMPASS), or - even better - 
if the gluonic 'background' can be eliminated or largely reduced by some suitable cuts. 
The latter can be partly achieved for the z integrated Z)-meson cross section d/\a /dx 
by introducing - as in Fig. 3 above - a lower cut-off z > 0.2 for the integration as can 
be inferred from Fig. 6. Contrary to the unpolarized case, the gluonic correction does 
not fully drop out here, since it is not oscillating in the region 0.2 < 2; < 1 (at least 
not for the chosen A^f). However, this observable nevertheless seems to be best suited 
for a determination of the strange sea, but some knowledge of /\g would be certainly 
desirable to disentangle its 'background' more precisely. It should be mentioned that the 
dependence on the factorization scale is rather weak, similar to our observations for the 
unpolarized case. 

Finally, let us study the sensitivity of CC charm production to the unknown spin- 
dependent strange density by comparing the results obtained for different, extreme choices 
of As. First of all it should be noted that the actual observable is the spin asymmetry 

for D meson production 

^ f dz dAa/dx/dz 
j^D ^ Jz^ '_ '__ (27) 

J dz da/dx/dz 



rather than the polarized cross sections for Zjam = and 0.2 shown in Fig. 6. A^ in (p7|) 
is simply related to a measurement of the counting rate asymmetry for parallel and anti- 
parallel alignment of the proton and lepton spins and does not require the determination of 
the absolute normalization. Furthermore, other experimental uncertainties conveniently 
drop out in the ratio (p7|). 



Fig. 7 shows our results for dAa/dx and A^ in longitudinally polarized e~p and e~^p 
collisions for z^m = and 0.2. Apart from the GRSV 'standard' set, we now adopt 
also the GRSV 'valence' set |2^, which, on the contrary, has a small positive strange 



density in the relevant x region whereas the gluon distribution is practically unchanged. 
While the cut z > 0.2 merely changes the size of the cross section and hardly effects the 
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asymmetry A^ for the 'standard' set, the influence of the cut on A^ is more pronounced 
for the 'valence' set as can be inferred from comparing Figs. 7(b) and (d). The oscillating 
behaviour in the 'valence' case stems from the interplay of Born and NLO contributions 
with different signs. Also shown in Fig. 7 is the expected statistical accuracy 6A^ for 
such a measurement at a polarized HERA 



5A- = j, , ^ ' (28) 

assuming an integrated luminosity of £ = 500 pb~^, 70% beam polarizations Pp and Pe, 
and where we have integrated over three bins in x. As can be seen the results for either e~p 
or e~^p collisions for the two different sets of parton densities can be distinguished within 
the error bars and hence some information on As can be extracted. However, there is a 
severe catch: In Fig. 7 a charm detection efficiency of 100% was assumed, i.e., e^jj- = 1, as 
was also used in previous LO calculations |^^, but which is completely unrealistic. With 



present-day values for e^rr at HERA of much less than 5% a measurement of As via CC 
D meson production is certainly impossible. However, until a polarized HERA could be 
realized in the future, some progress on charm detection might be possible. Furthermore, 
given the possibility that runs will be made with e~ and e~^ beams these results could be 
added to double the statistics. Of course, possible changes of sign in dAa have to taken 
into account by taking the absolute value, and it would be somewhat less clear how to 
extract As in such a case. 

4 Summary 

We have presented the complete NLO QCD framework for CC mediated inclusive charm 
and momentum z differential D meson production in DIS with unpolarized and longitu- 
dinally polarized beams. The calculated 0{as) coefficient functions refer to the MS fixed 
flavor scheme which fully takes into account the mass of the produced heavy (charm) 
quark. Our unpolarized results fully agree with previous calculations while the spin- 
dependent expressions are entirely new. Special emphasis was put on technical subtleties 
in the z differential case and due to the appearance of 75 in n dimensional regularization. 
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Exploiting our analytical results we have performed a detailed phenomenological anal- 
ysis of the prospects of determining the unpolarized and polarized strange density from 
CC D meson events at HERA. It was shown how to reduce the 'unwanted' contribution 
from the genuinely NLO gluon initiated subprocess considerably by imposing a lower cut- 
off on the D meson momentum fraction z for z integrated rates. Furthermore, it was 
argued that the sizeable NLO corrections observed for inclusive charm and momentum z 
differential D meson production can be both understood and even expected in the kine- 
matical domain of HERA due to the peculiarities arising from the mixture of weak and 
strong interactions in the case of CC charm production. 

It turned out that the rather small charm detection efficiency is the main limiting factor 
in pinning down the strange sea at HERA. Nevertheless, the unpolarized strange density 
was shown to be measurable with sufficient accuracy (a decent luminosity permitting), 
and HERA can hopefully contribute to a better understanding of the so far only weakly 
constrained strange distribution in the future. Unfortunately, in the polarized case, where 
much less is known about the flavor decomposition of the sea, a useful measurement of 
As from the CC D meson spin asymmetry cannot be performed without a significantly 
improved charm detection efficiency. 
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Appendix 



Here we list the expressions for the NLO (MS) CC coefficient functions AHf^F^ 4 ^ for heavy 
quark (charm) production appearing in Eqs. (pT]) and (^). The ( differential fermionic 
NLO (MS) coefficients AH^^^^^{^,(^ f^F^ ^) i^ (P^), as obtained from calculating the 
the subprocess W~^s -^ gc and the virtual corrections to W'^s —>■ c, coincide with the 



unpolarized functions -f^f=i 23(^5 C;/^F) '^) iii EH- They will be nevertheless also given 
here for completeness: 

Q^ + ml 



A/J|(e,C,/i^,A) = 5(1-0 {ap(J)(0 In : 



/^l 



4 

+ 3 



l-e + (l-01n 



;i-o^ 



+ -^-5(i-o<5(i-C) 



e(l-AO M-^ 



1 /1 + 3A ^^ 1 



2^ ^ + 2^(^ rln^ ^^ 



1-^ l-A^ 



4 + 



TT 



^ Az^^,i_0^i-^^ 



+ 2 



1-0. 



2 r 



+ 2e 



1-0+ (1-0 

l-Ae 



1-e 

i-(i-0 



1-e 

(l-AO^ 
l-Af 

i-e. 



i-(i-0 



i-e 



(Al) 



Ai7l(e,c,/^I,A) = Aff|(e,c,/i|,A) + ^|<5(i-05(i-o^A 

- 2(e(l-3A)[l-(l-C)^^] + (l-A) 



(A2) 

Ai7?(e,C,/^l,A) = AiJ|(e,C,/x|,A) +2^|(l-e)[l-(l-0\^]-(l-A0} (A3) 
where we have defined A = QV(Q' + ml) and Ka = }(1 - A) ln(l - A) and APig\^) = 



4 i±e 



3 V 1-5 



denotes the LO q —>■ q splitting function. The '+' and '©' distributions in 



1|)-(|A3|) are defined by 



\ ^'«' 



■,^ /«)-/(!) 



'0 (1-0+ io ""' l-( 
with C„,„ = (1 - \)(/(l - A{). 



Sniin 



:i - c)( 



' ,C ^'^,'-f' (A4) 

Smin 



1-C 



When integrated over (, the results given in 



reduce to the inclusive coeffi- 
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cients 

Ai7f (e, /x^, A) = /" dC AHfi^, C, f^% A) 

"^ Cmin 

appearing in (|2l|), where 
with 



and 



APW(0 In^^ + A/.^(e, 



/^F 



A) 



+ S; 



1 



2i 



;i - AO- 



+ ^3, 



0+ 



;i-Ao^ 



/i-? 






1-e 



21n(l-0-ln(l-A0 



(A5) 



(A6) 



(A7) 



(A^ 



The coefficients in (A7) for i = 3,4,1 are given in Table 1 and agree with the resuhs 
presented in p|, pi. 



Table 1. Coefficients for the expansion of A/i^ in 
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Bu 
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1-A^+e 
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^-e 



-1-^ 
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The C differential NLO (MS) gluonic coefficient functions AiJf (^, C,/i|,, A) for heavy 
quark (charm) production in (0), as calculated from the BGF subprocess W~^g —>■ cs, 
are given by 



AHl,A^,C,f^l,X) = <5(1-C)<|APW(0 

1 



+ 



1 



In 



Q^ + ml 



+ ln 



:i-o^ 



:i-c)e^c 



e(i-Ao, 

APi?(0 + A/.f(e,C,A) 



+ (1-0 



(A9) 
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where 



A/ii(e,C,A) 
A/iI(e,C,A) 



-e(l-A) 
1-A 



1 2 



^2 ^(1 - 2A) + - [2^(1 - A^) - (1 - A)] 



A/^?(e,C,A) = e(i-A) 



1 2 



+ 1 - 2Ae 



with APqg\^) = I [2^ — 1] denoting the LO q ^^ g sphtting function. The © distribution 
is aheady defined in ( |A4|) . When integrated over (, the results in ( |A9|) reduce to the 
inclusive coefficients in f^ll), i.e., 



Ai/f (e, /i|, A) = /" dc Ai/f (e, c, z^?., A) 

"^ Cmin 



(AlO) 



where 

Ai7^3,4(e,/i|,A) 

*~ 1 
with 



and 



Ai^J)(0 (tL, + In ^!±^ + In ^^3^) + A/^f (e, A) 

(All) 
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A/ii(e,A) = 2(l-A)L,e 

A;iKe,A) = 2(l-A)(l-0 + (l-A)L, [2e(l + A)-l] 

A/if(e,A) = (l-0[4-l/(l-AO]-2(l-A)L,e • 
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Figure Captions 

Fig. 1 Feynman diagrams contributing to the CC massive charm [p^ = m'^) production 
up to 0{as)'. Born term (a), real gluon emission (b), virtual corrections (c), and 
boson gluon fusion (d). The Cabibbo suppressed contributions are obtained by 
substituting all s hj d quarks. The relevant diagrams for s ^ c transitions are 
obtained by reversing all quark lines in (a)-(d). 

Fig. 2 z differential CC D meson production cross section e~p -^ DX in NLO obtained 
by integrating Eq. (|^) over the range Q^ = SepXy > 500 GeV^ with 0.01 < y < 0.9 



and ^/S^ = 300 GeV for HERA. The GRV-94 NLO parton densities H, m^ = 
1.5 GeV, fijp = Q^ +m^, and Ec = 0.06 in the charm fragmentation function ( |2^ ) are 



used. Also shown are the individual NLO gluon- and quark-initiated contributions, 
the latter including the virtual corrections, and the Born term obtained with NLO 
parton distributions. 

Fig. 3 The CC D meson production cross section e~p — > DX in NLO obtained from 
Eq. ([1^) as in Fig. 2 but now as a function of x and integrated over two ranges of 
z: < 2; < 1 and 0.2 < z < 1. Also shown is the Born contribution in each case, 
obtained with NLO parton distributions. 

Fig. 4 Comparison of the NLO z integrated D meson production cross section e~p -^ DX 



as in Fig. 3 using the GRV-94 [^ and MRST |1|| distributions for two ranges of z: 



{a)0< z<l and (b) 0.2 < z < 1. 

Fig. 5 As in Fig. 2 but now for longitudinally polarized e~p collisions according to 
Eq. ( |20D and using the GRSV 'standard' set of polarized parton distributions ||2^ . 

Fig. 6 As in Fig. 3 but now for longitudinally polarized e~p collisions according to 



Eq. (pOD and using the GRSV 'standard' set of polarized parton distributions |23 



Fig. 7 The z integrated polarized cross section [(a): < 2; < l,(c): 0.2 < z < 1] as in 
Fig. 6 and the corresponding measurable spin asymmetry A^ [(b): < 2 < l,(d): 
0.2 < z < 1] according to (|27|). Also shown for comparison are the results for 
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e+p -^ DX and the ones obtained using the GRSV 'valence' set of spin-dependent 
parton densities ||23[. The error bars denote the expected statistical accuracy 5A^ 



according to (pSD for three different x bins assuming C = 500 pb~^, Pe = Pp = 0.7, 
and eSr = 1. Two x bins are chosen at equal logarithmic distance in the range 
[xmin — 0.006, 0.1] and one for x > 0.1. 
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